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We have studied the shock-induced phase transition from fullerite to a dense amorphous carbon 
phase by tight-binding molecular dynamics. For increasing hydrostatic pressures P, the C6o-cages 
are found to polymerise at P < 10 GPa, to break at P ~ 40GPa and to slowly collapse further at 
P > 60 GPa. By contrast, in the presence of additional shear stresses, the cages are destroyed at 
much lower pressures (P<30 GPa). We explain this fact in terms of a continuum model, the snap- 
through instability of a spherical shell. Surprisingly, the relaxed high-density structures display no 
intermediate-range order. 

PACS numbers: 61.48.+C, 62.25.+g, 64.70.Nd, 02.70.Ns, 62.50.+p 



Rapid development in cluster science and nanotechnol- 
ogy has fostered the hope that materials exhibiting novel 
properties due to the existence of intermediate-range or- 
der (IRO) can be constructed from clusters Q or nan- 
otubes as building blocks. Fullerenes (Cgo molecules) 
represent one of the rare species of clusters available in 
macroscopic quantities Q • It is not surprising, therefore, 
that fullerene-assembled materials have been studied for 
at least a decade now 0. 

The extraordinary elastic properties of the molecu- 
lar cage of Ceo have inspired many to speculate about 
new carbon modifications with extraordinary mechani- 
cal performance [f|. Indeed, in early experiments very 
hard and stiff materials were synthesised (stable under 
ambient conditions) by compressing fullerite under high 
pressure 0. Subsequently, these phases have been char- 
acterised experimentally by determiningtheir structural, 
mechanical, and optical properties 0, H El CZJ EI 
EU Eil revealing a plethora of ordered (polymerised 
fullerenes) and disordered (amorphous) carbon phases 
(see for a review). 

Based on electron diffraction experiments, shock- 
compressed fullerite 0], has been conjectured to be a 
new form of amorphous diamond exhibiting IRO Q . Fur- 
thermore, it has been argued that the mechanical prop- 
erties arc determined by remnants of (partially) intact 
fullerene cages distinguishing these phases from tetrahc- 
drally coordinated amorphous carbon (ta-C) produced 
by ion-beam techniques. 

However, up to now, conjectures concerning the struc- 
tural properties arc based on indirect information; nei- 
ther the process of formation of such structures nor the 
nature of the proposed intermediate range order is under- 
stood microscopically. Pioneering molecular-dynamics 
(MD) simulations have shown that subjecting fullerite 
to high pressure and high temperatures may give rise 
to a dense amorphous phase yjj- However, the sample 
in 01 had to be compressed to very high densities (4.4 
g/cm 3 at T = 2500K) corresponding to pressures exceed- 



ing lOOGPa, far beyond any experimentally obtainable 
conditions. Empirically, the transition from fullerite to 
amorphous carbon occurs in the range of 10 — 30 GPa, 
depending on the speed of compression and on tempera- 
ture (room temperature in |6j and approximately 2000K 
in shock experiments HE3). How' may this contradic- 
tion be resolved? There are indications that shear stress 
may play a role in the destabilisation of the cages |[|. 
However, at present it is not known which microscopic 
processes cause the transition from fullerite to amorphous 
carbon. 

There is thus a pressing need for a theoretical under- 
standing of the exact microscopic mechanisms occurring 
during this phase transition. To which extent is a macro- 
scopic picture of a rapidly imploding spherical shell ap- 
propriate? Is shear important? How strong is the de- 
pendence of the final structural properties on the details 
of the compression process? Is there intermediate range 
order in the final dense carbon phase? 

This letter reports on quantum-mechanical computer 
simulations of shock-compressed fullerite, starting from 
fcc-crystallised and tetragonally polymerised fullerenes. 
Our findings provide profound insights into the above- 
mentioned issues and can be summarised as follows. 

Microscopic mechanisms and transition pressure. Hy- 
drostatic pressures transform fullerite into amorphous 
carbon by a succession of three steps: (i) polymerisation 
mainly by 2+2 ring-formation starting at P < 10 GPa, 
(ii) breaking of the Ceo cages starting at P ~ 40 GPa, 
and (Hi), at still higher pressures, collapse of the bro- 
ken Ceo cages via formation of chemical bonds across the 
cage. Remarkably, process (Hi) occurs very slowly, and 
(c.f. Ref. H3) very high hydrostatic pressures (exceeding 
60 GPa) are required for a complete collapse. 

Non-hydrostatic (external) stresses, by contrast, result 
in a rapid shear deformation of the unit cell: steps (ii) 
and (Hi) happen simultaneously and very rapidly. In 
other words, the Cgo molecules are ground - and not 
merely crushed to form amorphous carbon. Application 
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of shear significantly lowers the transition pressure. We 
explain this fact in terms of a continuum model (snap- 
through instability of a spherical shell under shear). 

Micro structure and intermediate-range order. We 
analyse the resulting structures after compression to 
35 GPa and subsequent relaxation to ambient conditions 
and find that the morphology depends on the amount of 
applied shear. In order to characterise the phase tran- 
sition and resulting morphologies it is not sufficient to 
merely specify pressure and temperature, shear is equally 
important. At low shear stresses the mechanical strength 
of the material is low and at the same time, intact Ceo 
cages are still present. At high shear stresses, by con- 
trast, no remnants of Gqq cages are observed. In this 
case, bulk moduli are found to be large, of the order of 
that of diamond. The radial distribution functions G(r) 
indicate that our fullerite-derived dense phase has lost its 
nanoscopic order and turned into conventional ta-C. We 
show that fine details in G(r) of shock compressed ful- 
lerite 0| are most likely effects of an incomplete exper- 
imental fc-space sampling and not a fingerprint of IRO. 

Methods. We simulate the MD of shock-compressed 
fulleritc within an fee and a tetragonal unit cell contain- 
ing four and two Ceo , respectively. Interatomic forces are 
determined within a non-orthogonal density-functional- 
based tight-binding scheme 0] . A Langevin thermostat 
is used to keep the temperature close to the 2000 K re- 
ported for shock conditions Q. Spontaneous shearing 
of the simulation cell is allowed for and external shear 
stresses were applied by using a Rahman-Parrinello 
barostat. After thermalisation of the sample, the exter- 
nal pressure is increased linearly from to 70 GPa within 
100 ps. Bond formation and breaking as well as sp 3 /sp 2 
hybridisation is determined based on a population anal- 
ysis of the bond orbitals, similar to the Mulliken charge 
analysis [2£j. In the remainder of this letter we explain 
how the results summarised above follow from our simu- 
lations (Figs. □ to EJ) • 

Hydrostatic conditions. Fullerite proves to be very sta- 
ble under isotropic compression, withstanding high pres- 
sures. The microscopic mechanism of collapse is appar- 
ent from Fig. displaying the evolution of inter- and 
intramolecular bonds and the representative fate of one 
of the cages. Starting from an fcc-ordcrcd crystal, the 
structure changes gradually through polymerisation of 
Cgo-molecules (see the rise in the number intermolecular 
bonds for P > 5 GPa, lower curve of Fig. la). Dou- 
ble bonds in neighbouring cages break and rchybridise 
with each other to form a mainly four-membered rings 
(2+2 cycloaddition (^IH^I). This process is completed 
for P > 26 GPa and the number of intermolecular bonds 
stays constant up to P c = 38 GPa (Fig. la). 

Upon increasing the pressure beyond P c = 38 GPa, 
the cages break (inset in Fig. la), indicated by a drop 
in the number of intramolecular bonds (upper curve in 
Fig. ^i) an d form new bonds to neighbouring fullcrenes 
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FIG. 1: Fullerite under increasing hydrostatic pressure, at 
2000 K. a Evolution of bonding in the absence of shear stress, 
starting from an fee crystal of fullerenes. Shown are the num- 
ber of bonds per C6o-molecule as a function of applied pres- 
sure. Intramolecular (red, upper curve) and intermolecular 
bonds (black, lower curve) are plotted separately. Insets show 
the shapes of a representative Ceo molecule as the pressure 
increases, in steps of 14 GPa, and the results of a continuum 
model (snap-through instability of a spherical shell) at pres- 
sures 5, 20, 30 and 40 GPa. The instability occurs between 
the third and fourth configurations, b Comparison of density 
evolution of the compressed Ceo crystal starting from the fee 
(black curve) and from a tetragonally polymerised phase (red 
curve). Insets show the initial and a intermediate configu- 
rations (at pressures and 35 GPa, respectively) for the fee 
(top row) and tetragonal (bottom row) phase. 

(sudden raise of the lower curve in Fig. [IJi) leading to 
a densification of the material (Fig. lb). At even higher 
pressures, the Cgo cages gradually collapse forming bonds 
inside the cages marked by a steady increase in the num- 
ber of intramolecular bonds in Fig. ^ f° r large values of 
P. When and how fast this process occurs depends on 
the local alignments of the cages. 

In order to study the influence of the initial configura- 
tion on the phase transition, a second simulation starting 
from a tetragonally polymerised state (lower inset in Fig. 
lb) has been performed. The resulting density-pressure 
curve (Fig. lb) is qualitatively unchanged at high pres- 
sures, as is the resulting structure (at pressures exceed- 
ing 60 GPa). In the intermediate regime, however, the 
structures are quite different (insets of Fig. QJd) , as are 
the density-pressure curves (Fig. ^d): in the tetragonal 
phase, the Ceo molecules cannot rotate freely (as opposed 
to those in fee phase). Intermediate configurations are 
thus seen to be very regular, with a high degree of poly- 
merisation. 
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FIG. 2: Fullerite compressed in the presence of external 
shear, a Evolution of the number of bonds per Ceo molecule 
(color coding as in Fig. la). Insets displays the morphology 
of a Ceo-molecule (top) as the pressure increases, in steps of 
6 GPa, and from a continuum model (bottom) . b Density as a 
function of pressure in the presence of external shear S = <J xy . 
Shown are four plots corresponding to four pathways in the 
5-P-plane as given in the inset: (i) black, (ii) red, (iii) blue, 
(iv) green. Locations of the collapse transition are shown as 
points in the S-P-plane. Also shown is the critical line (red) 
as determined from a continuum model. Finally, morpholo- 
gies of the C-atom configurations in a unit cell for case (iii) 
are also shown. 



Regardless of the initial conditions, unrealistically high 
pressures (exceeding 60 GPa) are required to entirely de- 
stroy the cages (resulting in a dense amorphous phase af- 
ter relaxation to ambient conditions). This is not consis- 
tent with the experimentally found transition pressures. 

Non-hydrostatic conditions. We have compressed the 
polymerised fee phase under non- hydrostatic conditions, 
along three paths - (ii), (Hi), and (iv) - in the S-P plane 
shown in the inset of Fig. 2b (S = a xy is the shear 
stress in the (001) plane). Path (i) corresponds to the 
hydrostatic case. Fig. 2a shows the numbers of intra- 
and intermolecular bonds per Ceo molecule, for case (Hi). 
Comparison with Fig. la and inspection of the pressure- 
density relationship in Fig. 2b confirms that shear plays 
an important role lowering the critical pressure and in- 
creasing the speed of the phase transition which hap- 
pens now within a few GPa at approximately 34, 31, and 
26 GPa for cases (ii), (Hi), and (iv). The critical shear 
values at the collapse are S c = 6.94, 11.8 and 12.5 and 
shown as points in the S'-P-plane in Fig. 2b. 



Load mode p [g/cm 3 ] B [GPa] % sp 3 



(i) 


2.1 


50 


18% 


(") 


2.3 


60 


20% 


(iii) 


3.0 


210 


69% 


(iv) 


3.1 


320 


79% 



TABLE I: Characterisation of fullerite phases after compres- 
sion to 35 GPa and relaxation to ambient conditions. Load 
conditions as in Fig. ^Jp. 

Fig. 2a also indicates that the microscopic mechanism 
is different from the zero-shear case described above: in 
the presence of shear the cages are ground (see the tran- 
sition from a spherical to an oblate cluster in the inset of 
Fig. 2a). After a rapid transition the number of intra- 
and intermolecular bonds remains constant (as opposed 
to Fig. [Tk). The cages are so strongly flattened that the 
number of intramolecular bonds cannot be increased by 
raising the pressure further. 

The final structure however (at pressures > 35 GPa) is 
similar to the structure obtained under hydrostatic con- 
ditions at pressures exceeding 60 GPa. 




FIG. 3: Structures of fullerite relaxed from from 
P=35 GPa,T = 2000 K to ambient conditions, a-b mor- 
phologies of the fee unit cells resulting from compression with 
two different shear values: (a) no shear and (b) shear ac- 
cording to path (iv) in the S-P-plane displayed in Fig. 2b. c 
radial distribution function G(r) averaged over 20 ps along 
a trajectory at ambient conditions: (i) sample after a hydro- 
static compression to P = 35 GPa, (ii) same compression but 
with high shear, (iii) reference sample without any interme- 
diate range order (see text), (iv) after hydrostatic compres- 
sion to P = 60 GPa, (v) experimental G(r) of shock-treated 
fullerite first , (vi) same curve as in ii, convolved with a peak- 
shape function sin(Qr)/(7rr) with Q = 12. OA -1 . 

A continuum model. To which extent can the phase 
transition in Figs. 1 and 2 be understood in terms of 
a continuum model? We have modeled a Cqq -molecule 
in fullerite as a spherical shell with diameter 0.71nm. 
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Its Young modulus (5.32 TPa) and effective thickness 
(0.08nm) arc determined from tight-binding simulations 
of a graphene sheet. Forces are applied at 12 nodes cor- 
responding to the fee symmetry of the fullerite crystal. 
The shell exhibits a snap-through instability at P c = 37.8 
GPa (see snapshots of the shell in Fig. la) to be com- 
pared with 38 GPa obtained above. It must be noted that 
the agreement is better than expected jl^ , and that the 
continuum model cannot explain the dynamics and mor- 
phologies at pressures beyond the collapse (see below). 

The decrease of the critical pressure with increasing 
shear stress can also be understood in terms of the contin- 
uum model (inset of Fig. 2a): shear destabilises a spher- 
ical shell under pressure and lowers the critical pressure. 
The continuum model predicts that critical shear stress 
and pressure are linearly related, the result is shown as 
a red line in the inset of Fig. 2b, in good agreement with 
the atomistic results for not too large values of shear 
stress. This shows that the onset of the instability of 
fullerite under shear stress can be described by elasticity 
theory. The ensuing collapse, however, must be modeled 
microscopically (in terms of breaking of bonds and rehy- 
bridisation). Finally, the continuum model fails at high 
shear stresses (c.f. path (iv) in the inset of Fig. 2b). 

Resulting morphologies. The question remains: what 
structures are obtained when the compressed fullcritcs 
are relaxed from - for instance - 35 GPa and 2000 K 
to ambient conditions? The morphology of the relaxed 
phase depends clearly on the amount of applied shear. 
Consider the relaxed structures shown in Fig. and b. 
They were obtained from samples compressed along the 
S-P-path (i) and (iv) (as shown in the inset of Fig. 2b), 
respectively. For no or low shear stress, densities, bulk 
moduli and sp 3 fractions are small (see load mode (i) and 
(ii) in Tab. 1), polymerisation is evident but the cages 
are still intact (Fig. Ek). At high shear (load mode (Hi) 
and (iv) in Tab. 1), by contrast, densities, bulk moduli 
and sp 3 fractions are much larger and no remnants of Ceo 
cages are found (Fig. 03). 

Does the radial distribution function reflect possible 
IRO in the relaxed samples? Curves i and ii in Fig. Et 
display G(r) for the two cases discussed above. In the 
zero-shear case, Ceo-specific distances are still present 
(the average nearest-neighbour bond length in graphite 
is 0.142 nm) thus showing clear evidence of IRO. In 
the other sample, by contrast, the radial distribution 
function is clearly different exhibiting the first peak at 
0.154 nm, as in diamond. Curve iii in Fig. 3c corre- 
sponds to a sample relaxed from P = 60 GPa, 5 = 
GPa to ambient conditions. Obviously this phase has a 
structure similar to ii. 

To test the hypothesis of a new form of amorphous di- 
amond, we have created reference sample exhibiting no 
IRO, obtained by isochorus heating of 240 carbon atoms 
to 10000 K while monitoring the system for sufficient 
chaotic motion followed by a relaxation to ambient condi- 



tions. It exhibits ta-C characteristics and approximately 
the same density and sp 3 /sp 2 hybridisation as the dens- 
est structure in Tab. 1 Surprisingly, the resulting G(r) 
(curve iv in Fig. |3t) is essentially identical to ii and con- 
sequently, the radial distribution function ii in Fig. 3b 
provides no fingerprint of IRO . 

The conditions of the simulations reported here are 
close to those in the shock-compression experiments re- 
ported in rcf. [lij . The third peak in the empiri- 
cally determined G(r) (curve v in Fig. Et) at approx- 
imately 0.312 nm has been accredited to the diame- 
ter of six-membercd rings in the Cgo molecules or to a 
slightly elongated third-neighbor distance for diamond. 
However, it should be kept in mind that experimental 
radial-distribution functions are generated from recipro- 
cal space scattering information which is limited by an in- 
evitable momentum cut-off Q. Thus, for the sake of com- 
parability, the simulated G(r) of our dense amorphous 
carbon (curve ii) was convoluted with an appropriate cut- 
off function reflecting the experimental Q=12 A -1 [lij . 
The resulting radial distribution function (curve vi in 
Fig. 3c) is very similar to that observed in [lfij |. We 
conclude that the third peak in v is unlikely to reflect 
IRO, it is likely to be an effect of the finite experimental 
resolution. 
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